Kinematic reconstruction of atmospheric neutrino events in a large water Cherenkov 

detector with proton identification 
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We report the development of a proton identification method for the Super-Kamiokande detector. 
This new tool is applied to the search for events with a single proton track, a high purity neutral 
current sample of interest for sterile neutrino searches. After selection using a neural network, we 
observe 38 events in the combined SK-I and SK-II data corresponding to 2285.1 days of exposure, 
with an estimated signal-to-background ratio of 1.6 to 1. Proton identification was also applied to a 
direct search for charged-current quasi-elastic (CCQE) events, obtaining a high precision sample of 
fully kinematically reconstructed atmospheric neutrinos, which has not been previously reported in 
water Cherenkov detectors. The CCQE fraction of this sample is 55%, and its neutrino (as opposed 
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to anti-neutrino) fraction is 91.7 ± 3%. We selected 78 ^-like and 47 e-like events in the SK-I and 
SK-II data set. With this data, a clear zenith angle distortion of the neutrino direction itself is 
reported in a sub-GeV sample of /i neutrinos where the lepton angular correlation to the incoming 
neutrino is weak. Our fit to Vr oscillations using the neutrino ^ distribution of the CCQE 

sample alone yields a wide acceptance region compatible with our previous results and excludes the 
no-oscillation hypothesis at 3 sigma. 

PACS numbers: 14.60.Pq, 14.60.St, 29.40. Ka 



I. INTRODUCTION 



detector. 



The Super-Kamiokande detector is well known for its 
discovery of atmospheric neutrino oscillations. It is a 
large, 50-kton, cylindrical water Cherenkov detector lo- 
cated under Mt. Ikenoyama near Kamioka, Japan [T]. 
It is separated in two concentric regions: a 2 m-thick 
outer veto (OD), instrumented with 1885 8 inch photo- 
multiplier tubes (PMTs) facing outward, used as a cosmic 
ray veto; and an inner detector (ID) instrumented with 
11,146 inward facing 20 inch PMTs. The Cherenkov light 
emitted by charged particles traveling in the water trav- 
els to the walls and the resulting ring-shaped patterns 
are then analyzed. 

Super-Kamiokande has very good particle identifica- 
tion (PID) capabilities, being able to separate muon 
tracks from electron tracks with better than 1.8% mis- 
identification [2]. Until now, no attempt has been made 
to extend its PID algorithm to other types of particles 
produced in neutrino interactions. In this article we 
demonstrate that Super-Kamiokande can also be used to 
search for proton tracks. Other authors have published 
estimates predicting what would be found in Super-K and 
other experiments, along with explorations of the physics 
consequences [3]. 

This improved PID technique will be described in sec- 
tion II B In section III we will describe atmospheric 
neutrino data selection and reconstruction at Super- 
Kamiokande. In section [IV| we will apply this new proton 
selection algorithm to a search for neutral-current (NC) 
elastic I'+p — > i^+p events. This sample has high neutral 
current purity and angular correlation with the incoming 
neutrino, and is potentially sensitive to i^steriie os- 

cillations [3]. 

In section |V] we will extend the technique to charged- 
current quasi-elastic (CCQE) neutrino interactions, re- 
constructing both the outgoing lepton and the proton. 
This allows full kinematic reconstruction of the atmo- 
spheric neutrino track, improving neutrino zenith angle, 
flight length and energy reconstruction, which is useful 
for oscillation studies. This is the first time such a re- 
construction has been performed in a water Cherenkov 
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II. A PROTON TRACK FITTER FOR 
SUPER-KAMIOKANDE 

A. Characteristics of proton tracks in 
Super-Kamiokande 

Proton tracks in water exhibit certain features that 
make it possible to identify them. 

The main parameter in Cherenkov light production is 
the velocity /3 of the particle, related to the Cherenkov 
cone opening angle by 



cos 6 
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where n is the refraction index of the medium. In wa- 
ter this leads to a threshold of /3 « 0.74. Because of 
the large proton mass this corresponds to a momen- 
tum threshold of 1.07 GeV/c. Due to the falling en- 
ergy spectrum of atmospheric neutrinos, we do not ex- 
pect to see many proton tracks above 3 GeV/c, corre- 
sponding to (3 < 0.96. Consequently for most of the pro- 
tons seen in Super-Kamiokande the half-opening angle of 
the Cherenkov cone will be relatively small, below about 
38° . This contrasts with electrons and muons which have 
lower mass, and so reach /? « 1 and Cherenkov angles of 
« 42° at much lower momenta than protons. 

Super-Kamiokande 's PID program uses the ring topol- 
ogy to separate muons and electrons: due to scattering 
and radiative processes creating electromagnetic showers, 
the outer edge of electron rings are blurred, while muons 
have rings with sharp outer edges [2]. Like muons, pro- 
tons will have rings with sharp outer edges. Since muons 
at these momenta lose their energy by ionization, their 
ranges are long (typically several meters) , leading to thick 
projected rings on the detector walls. However protons 
have a high probability of undergoing nuclear interactions 
with ^^O or H nuclei in the water. A large fraction of 
the time, neither the proton nor the secondary particles 
produced in these collisions are above Cherenkov thresh- 
old: in this case the proton's Cherenkov light emission is 
seen to stop suddenly, causing a short, thin ring to ap- 
pear on the detector wall. Figure [l] shows typical proton 
and muon event displays. 

Proton tracks have an additional characteristic: as 
shown in Fig. [2] the probability for a proton to produce 



FIG. 1: Event displays of a Monte-Carlo NC elastic event, with proton momentum of 1490 MeV/c (left), and a Monte-Carlo 
300 MeV/c muon. The proton stopped early, causing a thin ring pattern on the wall. The muon ring is thicker than most 
proton rings with similar opening angles. 



a TT^ or a charged secondary during a hadronic interac- 
tion in the water increases with momentum, reaching 
~ 50% at around 2 GeV/c based on our Monte-Carlo 
studies. In particular, if a 7r° is produced, at least 
one bright electron-like ring (from neutral pion decay 
gamma showers) will be seen in Super-Kamiokande, 
greatly reducing the chances of identifying the proton 
and accurately reconstructing the event. The fraction of 
protons that do not produce any visible secondary parti- 
cles goes from - 90% at 1.2 GeV/c to - 40% at 2 GeV/c. 
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FIG. 2: Probability of hadronic interactions in the water as 
a function of proton momentum. The clear region shows the 
fraction of protons that do not interact hadronically in the 
water. The cross-hatched region shows events whose inter- 
actions produce only sub-threshold secondaries. The region 
with a vertical hatch pattern corresponds to production of 
above-threshold charged secondaries, with no tt''. The hori- 
zontally hatched region shows the amount of 7r° production. 



B. Expected light pattern engine and fitter 

Particle identification is a hypothesis test, and is dealt 
with by calculating a ratio of maximum likelihoods. This 
is attempted after a previous fitter has already computed 
a vertex position, Cherenkov opening angle and candi- 
date ring direction. A light pattern engine generates the 
average Cherenkov pattern corresponding to the given 
vertex and track configuration inside the tank. Depend- 
ing on the hypothesis being tested, several parameters 
of the Cherenkov cone are then adjusted to maximize a 
pattern likelihood built from the observed and expected 
charges on individual phototubes. This procedure is car- 
ried out twice, using the two different particle types that 
are being studied to generate the light patterns. In this 
analysis, a special particle identification algorithm that 
takes into account proton-like features was developed. 
During the first stage, we assume the ring was created 
by a proton, and during the second stage it is assumed 
to be a muon. We will now briefiy describe the light pat- 
tern engine as well as the fit performed during the proton 
hypothesis test. 

1. Calculation of expected charge patterns 

Expected light patterns are produced using a program 
described in |4j. For proton identification we have ex- 
tended the program's capabilities. Using a GEANT4 
[5] Monte-Carlo simulator, we obtained tables of the 
Cherenkov photon density 

dr dp d cos 9 ' 

where r is the distance from the vertex, p is the particle 
momentum and 9 is the angle with respect to the particle 
track (cone axis). This table was built in "pure" water, 
i.e. all scattering and absorption effects were turned off; 
moreover, proton hadronic interactions and ^-ray produc- 
tion were turned off, but energy loss by ionization was 
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kept. Thirty batches of monochromatic protons at mo- 
menta from 1.1 to 4.0 GeV/c were used, with 50 MeV/c 
spacing until 1.5 GeV/c, 100 MeV/c spacing until 3500 
MeV/c and one batch at 4 GeV/c. The distance r was 
sampled every 50 cm. We used 400 bins in cos 9 from 1 
to 0.6 to sample the edge of the Cherenkov ring as finely 
as possible. 

For a given vertex, track direction and momentum p, 
the amount of Cherenkov light on a given phototube is 
calculated by linear interpolation from this table. The 
computed photon flux is azimuthally symmetric around 
the proton track. PMT solid angle and acceptance cor- 
rections, water attenuation and scattering effects are then 
applied to the light pattern. Since the photon flux is tab- 
ulated at flxed momenta, the Cherenkov opening angle 
of the cone is quantized in relatively large steps; conse- 
quently before any interpolation is done, each ^^'^ g 
distribution is reweighted so that its peak angle matches 
the Cherenkov angle at input momentum p to avoid any 
numerical problems. 

As explained earlier, one of the main features of 
protons is that their Cherenkov light emission can be 
stopped suddenly in the water. The pattern engine was 
modifled to produce light patterns for any path length: 
stopping the path length is simulated by masking the 
PMTs inside the inner edge of the cone. A smooth cubic 
spline is used to describe the effects of flnite PMT size 
for PMTs just on the edge of the lit region. 

A MINUIT-based flt [B] was built to flt the observed 
light pattern with an expected light pattern. The free pa- 
rameters in the fit are the proton momentum and proton 
track length. 

The expected light pattern then needs to be normal- 
ized. The relevant variable in Super-Kamiokande is Rtot, 
the amount of photo-electrons collected in a 70° cone 
around the track, correcting for water attenuation and 
scattering as well as PMT acceptance [7]. For electrons 
and muons, tables linking the particle momentum to Rtot 
are used for momentum determination throughout the re- 
construction. For protons, a one-to-one conversion func- 
tion does not exist since protons may stop abruptly in 
the water, yielding the same amount of light for differ- 
ent initial momenta. Using Super-Kamiokande's Monte- 
Carlo simulation with hadronic interactions turned off, 
conversion tables from proton momentum to Rt^t^ , the 
maximum achievable value, were calculated. The pro- 
ton average track length without hadronic interactions, 
Lmaxi was also obtained in this way. Then, by simulat- 
ing batches of protons at fixed path lengths, it was found 
that is only a function of ^ ^ and not of the initial 

momentum of the proton. Figure|3]shows this function as 
well as a fit. These calculations allow us to renormalize 
the expected light pattern to its expected Rtot{l,p) value 
for any given input length I and momentum p, ensuring 
that the maximum likelihood fit is well behaved. 

The fitter returns the maximum likelihood £proton, as 
well as the proton-like momentum and proton-like length 
of the event. In a separate stage, Lmuon, the likelihood 




FIG. 3: -Mwmc as a function of -j-^ — for Monte-Carlo batches 

"tot -i^max 

at different momenta and fixed path-lengths in the detec- 
tor. The red line shows the best fit to the functional form 
a(exp(-6;^^) - 1), with a = -1.29 and b = 1.5. The func- 
tion allows to calculate the expected amount of light for any 
proton pattern. 



that the pattern is muon-like, is calculated. 



C. Performance of the proton fitter 

The performance of this new algorithm has been 
tested using Monte-Carlo samples of muons, protons, 
and charged pions. The samples were produced using 
GEANT3 and include full fluctuations of all processes 
and other detector effects. The momentum resolution for 
protons varies from about 3% below 1.5 GeV/c, to about 
10% at 2 GeV/c, reaching 30% at 2.5 GeV/c. The mean 
resolution of the momentum fit over the NC elastic pro- 
ton spectrum (for the atmospheric neutrino spectrum) is 
about 7%, with a small positive bias in the momentum 
determination of 3%. The resolution gets worse as the 
momentum increases because it is mostly controlled by 
the opening angle measurement of the Cherenkov ring. 
As the angle increases, small variations in the angle mea- 
surement result in larger changes in the fitted momen- 
tum. Figure [4] summarizes the results of momentum re- 
construction. 

Figure |5] shows the distributions of the fitted path 
lengths for mono-energetic simulated protons for two dif- 
ferent momenta. At low energy, the distribution has a 
sharp peak at higher values, corresponding to protons 
that traveled their full path-length (no interaction in the 
water), and a region at lower fitted track lengths corre- 
sponding to proton tracks that stopped early. At higher 
input momenta all protons traveled less than their max- 
imum path length and the distribution falls off. This 
observation of hadronic interactions in the water can be 
summarized as the average ratio — as a function of the 
input value of (ij. Figure |6] shows that the ratio clearly 
decreases with input proton momentum, showing the in- 
crease of the hadronic cross-section with momentum. 
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FIG. 4: Bias (left) and RMS (right) of proton momentum 
reconstruction, using monochromatic Monte-Carlo. 
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FIG. 5: Fitted track length for Monte-Carlo protons at 1.3 
and 2.0 GeV/c. At 2.0 GeV/c the protons never reach their 
maximum path length. 



III. THE SUPER-KAMIOKANDE 
ATMOSPHERIC DATA SET 

A. Fully-contained atmospheric data 

The data used in this paper were collected at Super- 
Kami okande during two separate run periods. The first 
one, labeled SK-I, ran from April 1996 to July 2001, with 
40% photocathode coverage (11,146 ID PMTs). The sec- 
ond one, called SK-II, ran from September 2002 to Oc- 
tober 2005, with 19% photocathode coverage (5,183 ID 
PMTs) , as a result of partial reconstruction following an 
accident that occurred in Nov 2001. In SK-II the PMTs 
were encased in a protective fiberglass shell with a trans- 



it 2 3 

p y = p /Mc 

FIG. 6: Average fraction of the maximum path length in 
water traveled in water as a function of proton momentum 
(estimated with NC elastic Monte-Carlo events). The effect 
of hadronic interactions is clear. 



parent acrylic window to protect against the risk of tube 
implosion. The two data sets respectively correspond to 
1489.2 and 798.6 integrated days of livetime. 

Super-Kamiokande data first undergo a reduction pro- 
cess, which selects atmospheric neutrino events and splits 
them into separate samples and eliminates backgrounds. 
In this paper, the only relevant sample is the fully- 
contained (FC) data set, since protons have relatively 
short tracks in the water and do not escape the tank. 
The purpose of FC reduction is to remove events entering 
the detector - mainly downward-going cosmic ray muons 
- as well as events caused by flashing PMTs, and events 
with tracks escaping the inner detector. The optical sep- 
aration between the inner detector (ID) and the outer 
detector (OD) is used during this process. Details of the 
reduction process can be found in [2]. 



B. Event reconstruction 

Fully-contained data is then fitted by a series of pro- 
grams which find the vertex and particle tracks. Details 
about the reconstruction can be found in [2j |7]. The 
initial step is to obtain a starting vertex and track direc- 
tion using timing information only. PMTs on the edge 
of the brightest Cherenkov ring are then identified. The 
method for tagging edge PMT, and Cherenkov opening 
angle measurement, is different in SK-I and SK-II. In SK- 
I, the distribution Q{0) of observed charges as a function 
of opening angle is constructed. The edge of the ring is 



taken to be the first zero of the second derivative 
occurring after the maximum of Q. In SK-II, with half 
the photo-coverage, the sampling of Q{0) was found to 
be too low for an accurate determination of the deriva- 
tives; therefore, a method based on light patterns is used. 
The expected light pattern of a Cherenkov ring (e or /i) 
is fit to the observed pattern, allowing its opening angle 
to vary; the edge of the ring is found at the best fit, and 
edge PMTs are then tagged. 

Once edge PMTs are tagged, they are used to refine 
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the Cherenkov ring direction. An estimated track length 
is calculated, and the vertex position is changed by max- 
imizing an estimator based on PMT timing and charge, 
treating PMTs inside and outside of the Cherenkov edge 
differently. This procedure is iterated until a stable ver- 
tex and direction are found. The next step is to search 
for other rings in the event. 

A Hough-transform method [5] is used to look for possi- 
ble ring candidates. A pattern likelihood function is built 
to determine whether adding each Hough seed matches 
the observed charge pattern better. Other estimators in- 
volving the charge density in the candidate ring are also 
considered. If a better match is found by adding a ring, 
the seed is added to the list of tracks. This procedure is 
iterated until a maximum of 5 rings is found. 

The next step is particle identification (PID): each ob- 
served ring is compared to a muon pattern and an elec- 
tron pattern using a likelihood method, and is then as- 
signed a particle type, either e-like or mu-like. After 
PID, single ring events are processed by a refined track 
fitter that uses PID information to improve the vertex 
and track fit. 

An important variable is the visible energy, which is the 
energy of an electron that would cause the same amount 
of light in the tank. It is estimated at the end of the 
reconstruction process. 

The performances of PID and the refined vertex fitter 
are not optimal for proton tracks, since by design the 
particle is assumed to be a lepton in these algorithms. In 
the search for NC elastic events described in the present 
section we do not use the results of these two algorithms. 
However, for the analysis described in section |V] we do 
apply them. 

C. Monte-Carlo simulation 

Monte-Carlo simulation of atmospheric neutrino 
events in Super-Kamiokande relies on three successive 
steps: 

1. Atmospheric neutrino flux determination: the 
model by Honda et al.[9 is used in this paper. 

2. Simulation of neutrino interactions in the water: 
Two different programs were used: the NEUT sim- 
ulation [To] and the NUANCE (version 3) [11] sim- 
ulation. Details about NEUT are also available 
in [5]. The versions used in this paper includes 
rescattering of the outgoing nucleons, which is rele- 
vant for estimating the number of protons expected 
above Cherenkov threshold. Both NEUT and NU- 
ANCE assume the axial- vector form factor to have 
a dipole shape, with the axial mass set to 1.21 GeV 
(for CCQE and CC single pion interactions). 

3. Tracking of interaction products in Super- 
Kamiokande : A complete simulation of Super- 
Kamiokande using the GEANT3 [12 package is 



used. The Cherenkov light production and propa- 
gation models have been tuned to calibration data 
[T]. In addition, careful studies of the reflection 
of light on the acrylic surfaces, PMT glass and 
plastic lining of the tank have been carried out. 
The output of the GEANT3 simulation is in the 
same format as the data, and undergoes the same 
reduction and fitting process. The events are 
reweighted for livetime, solar wind activity and for 
— !■ neutrino oscillations. 



IV. SEARCH FOR NC ELASTIC u + p^u + p 
EVENTS IN SUPER-KAMIOKANDE 
ATMOSPHERIC DATA 

The first application of the proton fitter described 
above is to search for single-ring proton events caused 
by the neutral current elastic interaction v + p —> v + p. 



A. Single proton selection 

We initially select a sample of fully-contained atmo- 
spheric events, with their vertices (calculated by the ini- 
tial fitter) more than 2 m away from the wall, which 
defines a 22.5 kton fiducial volume. We also require that 
there be only one ring. The vertex resolution - defined as 
the radius of the sphere centered on the true vertex that 
contains 68.3% of reconstructed vertices - for single pro- 
ton tracks from simulated atmospheric NC elastic events 
is estimated to be ^ 85 cm, with an angular resolution 
of ~ 2.8° on the reconstructed direction. The Cherenkov 
angle is measured with an accuracy of ~ 1.6°, with a bias 
of -1.4°. 

We estimate that 7.0 visible NC elastic interactions 
occur every year in the fiducial volume. These are events 
for which proton NC elastic collisions produced a single 
outgoing proton above threshold which in turn produced 
no visible secondary in hadronic collisions in the water. 

However, other interaction modes create similar events 
with a single visible proton, especially neutral current 
single-pion (NCItt) events where the pion is absorbed; 
and NC elastic collisions on neutrons where the neutron 
produces a proton above Cherenkov threshold through 
hadronic collisions. Using Monte-Carlo information, we 
determined that 5.9 such events are produced per year 
in the fiducial volume, for a total of 12.9 single visible 
protons per year. Of those, 11.4 (88%) events are fitted 
as one ring events with their vertex inside the fiducial 
volume by our reconstruction. In the remainder of this 
section we explain how to separate these fitted events 
from the large backgrounds using the proton identifica- 
tion tool described in section IIIBI 
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1. Removal of low energy backgrounds 

In atmospheric neutrino analyses, data with less than 
30 MeV of visible energy are discarded; they are not used 
for oscillation searches. However in the present analysis 
such a reduction cut is too stringent and would remove 
about 20% of visible protons. At visible energies below 30 
MeV, the FC atmospheric neutrino sample does contain 
neutrino interactions, but also low energy backgrounds, 
especially spallation events. These are radioactive (3 de- 
cays of fragments left over from the collisions of cosmic 
ray muons with ^^O nuclei in the water. We use the spal- 
lation removal technique used in solar neutrino searches, 
which is described in detail in |13j . Spallation events are 
correlated in time and space with one of the previous 
cosmic ray muon tracks; there is also a correlation with 
muon energy. 

As in the solar analysis, a likelihood function based 
on these three parameters was built. The vertex position 
and event time of each FC event is checked against all the 
cosmic ray muons occurring during the previous 100 s, 
and the highest value of the likelihood is saved. For SK- 
I, two different likelihood functions are used, depending 
on the success or failure of the cosmic ray muon track fit. 
For SK-II both likelihoods were merged into one. Cuts 
on the likelihood value were chosen to remove the events 
with a high correlation with a preceding muon. This cut 
was only applied if E^s < 50 MeV. Using real muon data 
and uniformly distributed times and vertices, we deter- 
mined that these cuts incur a 8.2% (9.1%) inefficiency on 
non-spallation atmospheric neutrino events in SK-I (SK- 
II) below 50 MeV. We also removed 2.7 days of livetime 
from the SK-I sample when performing the spallation cut 
because of missing cosmic ray muon data. The value of 
the cuts and induced inefficiencies differ from those in 
[13j because we are looking at atmospheric events after 
fully-contained reduction and vertex reconstruction, and 
not low energy data after solar neutrino reduction. 

In addition, an extra cut was developed to remove some 
remaining low energy events with a sparse ring, i.e. an 
incomplete Cherenkov ring pattern. Using the fitted ver- 
tex and direction, an 6°-wide annulus was constructed 
around the Cherenkov ring, and binned in the azimuthal 
direction in 36 (/)-hms. A sparse ring will have empty 
0-bins, and the average charge per bin will be low. A 
two-dimensional cut on events with an average charge of 
less 1 (0.5) pe and more than 15 (22) empty 0-bins in 
SK-I (SK-II) was applied. Studies show that this cut im- 
proves signal to background significance at later stages. 
Good agreement between data and MC is reached once 
the low energy backgrounds are removed, as can be seen 
in Fig.[7l 



2. Proton track selection 

After the spallation and sparse ring removal cuts, three 
selection criteria are applied on the sample to select pro- 
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FIG. 7: Visible energy distribution for data and MC single- 
ring FC events for SK-I and SK-II is shown with no spallation 
cuts applied (crosses), after the spallation cut (open circles), 
and after the sparse ring cut (triangles). The corresponding 
MC distributions are shown (full line: no cuts, dash-dotted 
line: only the spallation inefficiency is applied, thick dashed 
line: the sparse ring cut and spallation inefficiency are ap- 
plied). Spallation events are present only in the data. 



ton events: 

1. Events with more than 200 MeV of visible energy 
are discarded. This removes the bulk of the events 
used for the usual atmospheric analyses, caused by 
high energy muon or electron tracks. 

2. Events with a Cherenkov opening angle greater 
than 37° (35° in SK-II) are rejected. This effec- 
tively removes high energy muons, as well as show- 
ers caused by electrons or 7 rays from tt" decays. 

3. Finally, the pattern likelihood difference calculated 
by the proton fitter described in section |ll] viz. 

AlogL = log(ipi.oton) - log(imuon) 

is required to be positive. This selects events with 
a proton- like pattern. 

The distribution of the variables are shown in Fig. [8] 
The efficiencies of all cuts are summarized in tables |l] and 
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[nj In these tables, events labeled as signal are events with 
a single visible track caused by a proton above Cherenkov 
threshold; this is comprised of ^ 55% NC elastic events 
on protons, the rest being mostly neutral current pion 
production collisions where the pion was absorbed. These 
interactions are an irreducible background to a search 
for proton NC elastic events, but we include them in the 
signal since this analysis technique is designed to maxi- 
mize the efficiency of proton track selection while reduc- 
ing contamination from muon and pion tracks. At this 
stage, there are ~ 25.6 (~ 10.4) expected signal events 
in SK-I (SK-II). 



2. Neural network response 

Figure [lO] shows the response of the neural networks, 
labelled NN in what follows. Using the Monte-Carlo, 
we have studied the composition of the signal and back- 
ground as a function of neural network response NN. 
Figure [TT] shows the results. It should be noted that 
the neutral current content of the sample before any NN 
cut is 66% (45%) for SK-I (SK-II); when keeping only 
events with NN > 0, this fraction reaches 85% (76%). 
The NN > cut was chosen because it maximizes the 
significance signal/i/signal -I- background. 



B. Neural network selection 

At this stage, the signal-to-noise ratio has been de- 
creased from ^ 1/200 to ~ 1/3. It is possible to im- 
prove this using multivariate analysis techniques. Here 
we chose to use a neural network. We used the TMVA 
library for this calculation. 



1. Neural network input variables and training 

As inputs to the network, seven variables were used: 
the Cherenkov opening angle, the /i-like momentum of 
the event, the fitted proton-like length, the fitted proton- 
like momentum, the pattern likelihood difference, the 
number of decay electrons in the event, and the normal- 
ized length. The normalized proton length L„ is defined 

by 



where L and P are the measured proton-like track length 
and momentum, and LmaxiP) is the maximum path- 
length in water for a proton of momentum P. Figure [9] 
shows the distributions of these variables for data and 
Monte-Carlo. Signal events are shown as hatched his- 
tograms. 

A large sample of NC elastic events on protons was 
prepared to be used as signal for training purposes. Sev- 
eral samples of low momentum monochromatic charged 
pions (from 180 to 1000 MeV/c) and muons (from 140 to 
500 MeV/c) were used as background training samples. 
All reconstruction algorithms and previous selection cuts 
were applied to those events, which were then fed into the 
neural network for training. Checks were conducted to 
make sure no overtraining occurred. The best architec- 
ture was found to be a single hidden layer with 6 neurons. 

Since the distributions of the discriminating variables 
are different for SK-I and SK-II, two separate networks 
with identical architectures were used. 



C. Results from the high purity neutral-current 
proton enriched sample 

After the NN > cut we expect 34.3 events in SK-I 
data and 15.3 in SK-II data; out of those, 22.2 and 8.5 
respectively are events with a single track caused by a 
visible proton; 13.1 and 5.0 (respectively) are neutral cur- 
rent elastic collisions on protons according to the NEUT 
simulation. We observe 27 events in SK-I and 11 events 
in SK-II. 

We calculated the of the data neural output distri- 
butions to Monte-Carlo with and without single visible 
proton events, (using 6 bins from -1 to 1 to increase the 
bin contents for statistical purposes). Using the Monte- 
Carlo including all events the is 9.3 for 6 bins (prob- 
ability 15.7%); when all visible protons are excluded the 
value is 15.8 (probability 1.5%). The data is clearly 
compatible with the observation of single visible protons. 

We have determined that in this final sample the 
proton-like particle has a mean angle of 38° (resp. 41° 
for SK-II) with the incoming neutrino track, which means 
that the proton still has some correlation with the orig- 
inal neutrino direction and allows zenith angle studies. 
In order to cancel out any fiux normalization uncertain- 
ties we calculate the up-down asymmetry where 
up-going events are those with — cos^zonith < —0.2 and 
down-going events have — cos^zenith > 0.2. The results 
are summarized in table |IV| and Fig. |12| and are com- 
patible with no sterile oscillations. Clearly the statistics 
of our sample are too weak to extract meaningful infor- 
mation on Z/^ Vsterile OSCillation. Full l^sterile 

oscillations (already ruled out by [15]) would result in 
an asymmetry of « —15% which is comparable with our 
data's present statistical error. Super-Kamiokande has 
also already constrained the admixture of sterile neutri- 
nos to be less than 26% at 90% CL[T6] using the so-called 
2+2 model parametrized in |17j . The amount of sterile 
admixture is not measurable with the single visible pro- 
ton sample alone. Although such a study is beyond the 
scope of this paper, the single visible proton sample could 
be merged with other samples sensitive to sterile neutri- 
nos and included in a larger fit. 
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FIG. 8: Distribution of the visible energy, Cherenkov opening angle, and proton pattern likelihood for single ring events after 
spallation and sparse ring removal, for SK-I (top row) and SK-II (bottom row). The full and dashed line show Monte-Carlo 
expectation from NEUT and NUANCE, with oscillation reweighting, but no correction for the absolute flux normalization. 
The hatched histogram shows the fraction of signal events according to NEUT. 



Super-Kamiokande-I 


Data 


Total MC 


Signal MC 


Total MC 


Signal MC 




NEUT 


NEUT 


NUANCE 


NUANCE 



EC, FV, single-ring, spallation removed 8946 (100%) 8138.1(100%) 45.1 (100%) 8031.5 (100%) 41.2 (100%) 

Sparse ring removal cut 8509 (95.1%) 7729.7 (95.0%) 31.7 (70.4%) 7673.4(95.5%) 29.3 (71.1%) 

< 200 MeV 2101 (23.5%) 1894.2 (23.3%) 29.7 (65.9%) 1843.5 (23.0%) 27.9 (67.7%) 

Cone opening angle < 37° 1161 (13.0%) 1020.0 (12.5%) 28.9 (64.2%) 1009.4(12.6%) 26.6 (64.5%) 

Pattern ID estimator cut 74 (0.83%) 68.8 (0.85%) 25.6 (56.8%) 65.8 (0.82%) 22.7 (55.0%) 



TABLE I: Summary of SK-I NC elastic selection. In this table the Monte-Carlo has been reweighted according to live time, 
solar activity for SK-I, as well as neutrino oscillations assuming Am^ — 2.5 x 10~^ eV^ and sin^ 29 = 1. Both NEUT and 
NUANCE simulations return similar results. 



V. SEARCH FOR CHARGED-CURRENT 
QUASI-ELASTIC EVENTS 

Another study made possible with proton identifi- 
cation is the search for charged-current quasi-elastic 
(CCQE) events with both the lepton and the proton 
above Cherenkov threshold. The reaction is + n 
p + lepton. Having a proton in the final state can 
only happen for neutrinos and not anti-neutrinos, pro- 
viding a possibility for neutrino tagging, which has so 
far not been done in Super-Kamiokande. Knowledge 
of both outgoing tracks allows kinematic reconstruction 



of the incoming atmospheric neutrino, which has never 
been reported in a large water Cherenkov detector: in 
the zenith-angle Super-Kamiokande analyses with atmo- 
spheric fully-contained data [5], the lepton direction in 
the data and Monte-Carlo are compared directly. From 
our Monte-Carlo simulations, we expect about 150 (75) 
events in the SK-I (SK-II) data sample with both tracks 
above Cherenkov threshold, assuming 2-fiavor oscilla- 
tions with (Ato2 ^ 2.5 X IQ-^ eV^sin2 26' I). 

We have focused on events with either one or two rings 
reconstructed by the standard ring finder. Each class 
requires a different analysis technique and will be studied 
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Super-Kamiokande-II 


Data 


Total MC 


Signal MC 


Total MC 


Signal MC 




NEUT 


NEUT 


NUANCE 


NUANCE 



FC, FV, single-ring,spallation removed 4700 (100%) 4190.1 (100%) 18.9 (100%) 4174.2 (100%) 17.9 (100%) 
Sparse ring removal cut 4464 (95.0%) 4004.6 (95.6%) 13.9 (73.3%) 3992.0(95.6%) 12.6 (70.4%) 
E^is < 200 MeV 1169 (24.9%) 969.4 (23.3%) 12.8 (67.5%) 953.0 (22.8%) 11.8 (66.1%) 
Cone opening angle < 35° 261 (5.55%) 217.1 (5.18%) 11.6 (61.1%) 216.7 (5.2%) 10.6 (59.5%) 
Pattern ID estimator cut 37 (0.79%) 39.9 (0.95%) 10.4 (54.9%) 39.6 (0.95%) 9.7 (54.5%) 



TABLE II: Summary of SK-II NC elastic selection. In this table the Monte-Carlo has been reweighted according to live time, 
solar activity for SK-II, as well as neutrino oscillations assuming Am^ = 2.5 x 10^'^ eV^ and sin^ 26 = 1. The cut on the cone 
opening angle is tighter than in SK-I, since the fitter relies on a different ring edge finding method. Both NEUT and NUANCE 
simulations yield similar results. 



Run period 


Data 


Expected signal 


Expected background 


SK-I 


27 


22.1 


12.2 


SK-II 


11 


8.5 


6.8 



TABLE III: Summary of neural network selection. 



Up Down ± Stat 

SK-I Data 11 11 ± 0.21 

SK-I NEUT MC 13.76 14.44 -0.024 ± 0.026 

SK-I NUANCE MC 12.61 11.82 0.032 ± 0.038 



SK-II Data 


2 


5 


-0.43 


± 


0.34 


SK-II NEUT MC 


5.88 


6.38 


-0.041 


± 


0.033 


SK-II NUANCE MC 


5.44 


6.23 


-0.068 


± 


0.040 


Combined Data 


13 


16 


-0.10 


± 


0.19 


Combined NEUT MC 


19.63 


20.81 


-0.029 


± 


0.021 


Combined NUANCE MC 


18.05 


18.05 


0.00 ±0 


.029 



TABLE IV: Up and down-going events with corresponding 
asymmetries. All NC elastic selection cuts have been applied, 
including the A^A^ > cut. In this table the Monte-Carlo has 
been reweighted according to live time, solar wind activity 
for SK-II, as well as neutrino oscillations assuming Am^ = 
2.5 X IQ-^ eV^ and sin 29 = I. The error bars for simulated 
events correspond to Monte-Carlo statistics. 

separately. Figure [13] shows two typical CCQE events, 
one for each class. 



A. Events with two fitted rings 

The first case to consider is that for which both rings 
were identified by the standard ring finder described ear- 
lier. When this happens, in more than 97% of the cases 
the proton is found as the second ring, so we make that 
assumption in this analysis. The vertex, ring directions 
and lepton momentum are known from the reconstruc- 
tion algorithms. We have modified the light pattern en- 
gine to superimpose a lepton and a proton pattern; the 
lepton ring is kept fixed, its parameters being set to the 
fitted information. The likelihood of this combined pat- 



tern is then maximized using the same MINUIT-based 
fitting, varying the proton momentum and track length. 
The normalization of the light patterns is done using the 
same procedure as for single ring protons, simply adding 
the lepton pattern to it. 

Then, a second likelihood corresponding to a lepton 
and a muon light pattern is calculated. The discrim- 
inating variables are the fitted proton-like momentum 
and track length, as well as the log-likelihood difference 

log Z/icpton+proton lo§ -^Icptoii+muon ■ 

CCQE selection relies on the following cuts: 

1. We first select only vertices inside the fiducial vol- 
ume, with 2 fully-contained rings, and 100 < i?vis < 
4000 MeV (Note that for the CCQE search the en- 
ergies are high enough that no spallation cut is nec- 
essary) . 

2. Our next cut aims at removing background coming 
from pion production, where the second ring was in 
fact caused by a charged pion. In this case, the fit- 
ted proton-like momentum of the track is typically 
about 1.6 GeV/c, because the charged pions have 
typical momenta between 200 and 300 MeV, and 
thus a ring opening angle of roughly 30°. However, 
the fitted path length is very short (about 20 cm) 
because of pion interactions. In order to remove the 
charged pions, a linear cut on the proton-like mo- 
mentum and track lengths A < P < B x L{m) + C, 
with A = C =1.1 GeV/c and B = 7/6 GeV/c/m. 
This selection can be seen in Fig. [14] 

3. We then require that the opening angle of the sec- 
ond ring be less than 34°, to select particles with 
low p. 

4. To select proton tracks a cut on the pat- 
tern likelihood difference log(Liepton+proton) - 
log (iiopton-i- muon) > is applied. 
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FIG. 9: Discriminating variables used as input for the neural network. The expectations from NEUT (full line) and NUANCE 
(dashed line) are shown. The hatched histogram shows the fraction of signal events according to NEUT. 



5. At this stage the sample still contains non-CCQE 
events, mostly corresponding to mis-identified pi- 
ons. To improve CCQE selection we use event kine- 
matics. Let V be the 4- vector 

V^Pp + Pl-Pn, 

where Pp, Pi, and P„ are the 4- momenta of the 
proton, lepton, and target neutron. For a true 
CCQE event, V must be equal to Pi,, the neu- 
trino 4-momentum, and consequently the Lorentz 



invariant quantity V'^ must be ~ eV^/c"*. 
On the other hand, for a non-CCQE event, V is 
not equal to the neutrino 4-momentum and its 
square will be non-zero. In practice, nuclear effects 
such as Fermi motion of the neutron or scattering 
of the outgoing nucleon, and detector effects will 
smear the resolution on V'^, but it is still a pow- 
erful tool, as can be seen in Figs. [T4| and [TSl with 
CCQE events peaked around 0. We apply a cut 
-0.75 < y2 < 1.5 GeV/c^ for 2-ring events. 
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FIG. 10: Neural network response for data and MC (full line: NEUT, dashed line: NUANCE). Left: SK-I, Middle: SK-II, 
Right: combined. The MC were oscillated with (Am^ = 2.5 x 10"^ eV^, sin^ 26 = 1.0). 
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FIG. 11: Top: SK-I neural network response for signal and background, separated according to the main interaction modes. 
For signal events, the main mode is NC elastic on proton, followed by NC single pion production where the pion is absorbed. 
The rest is mostly NC elastic on neutrons where the neutron produced a proton. Background mostly comes from low energy 
CCQE muon events, and pion tracks from single pion interactions. 



Figures 14 and 15 show the distributions of these vari- 
ables for data and Monte-Carlo for SK-I and SK-II, show- 
ing good agreement. The results are shown in tables [V| 
and I VI I The overall signal efficiency is 42-45%, with an 
expected signal-to-background ratio of 1.4 to 1. The high 
signal efficiency can be explained by the fact that those 
protons have high enough momenta to have been found 
by the ring fitter, and are well above Cherenkov thresh- 
old. Monte-Carlo studies of the selected events show that 



of the 38 (19) events in the final SK-I (SK-II) sample, 33 
(16) do have their second ring caused by a proton (the 
rest mostly come from gammas and charged pions), and 
22 (11) of those are CCQE signal events. The non-CCQE 
proton rings mostly come from CC pion events, and are 
an irreducible background to the CCQE search. It is 
worth noting that at this stage the sample contains 93% 
neutrinos and 7% anti-neutrinos. 
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FIG. 12: Up-Down asymmetries with statistical errors. In each plot the present measurement is the left-most point, labelled 
"data". The middle-points, labelled Neut and Nuance, show the MC expectation including 2-flavor i/^ —> oscillation, with 
MC statistical error only. The right-most point labelled "sterile" shows the expectation of the Neut MC assuming oscillations 
are only Vp. VsterUe- The error bars for simulated events correspond to Monte-Carlo statistics. 








FIG. 13: Event displays of Monte-Carlo CCQE events, where both the muon (large ring) and the proton (smaller ring) are 
visible. In the left figure, both tracks were found by the standard fitter (the thin lines show the fit results). In the figure on 
the right only the lepton (a muon in this example) was fitted; the proton is visible as a smaller ring next to the muon. 



SK-I Data Total MC Signal MC 

FC, FV, two-ring, 0.1 < E^is < 4 GeV 1876 (100%) 1773 (100%) 50.12 (100%) 
l.K P < 7/6L(m) -I- 1.1 GeV/c 829 (44.19%) 839.4(47.34%) 36.12(72.1%) 
Angle of 2"^^ cone < 34° 135 (7.2%) 124.1(7.0%) 27.2 (54.3%) 

Pattern ID cut 86 (4.58%) 82.88 (4.67%) 25.26 (50.4%) 

-0.75 < < 1.5 GeV^ 44 (2.35%) 38.59 (2.18%) 22.71 (45.31%) 

TABLE V: Summary of SK-I CCQE selection for two-ring events. In this table the Monte-Carlo has been reweighted according 
to live time, solar wind activity for SK-I, as well as neutrino oscillations assuming Am^ = 2.5 x 10"'^ eV^ and sin^ 20 = 1. 



B. Events with a single fitted ring after standard 
algorithms 



Monte-Carlo studies show that there are potentially 
w 150 CCQE events in the SK-I and II data sets with a 
proton above threshold where the proton was missed by 



the standard ring finder. Such events only have a sin- 
gle fitted ring corresponding to the brighter lepton track. 
The proton ring is weak, but visible by eye-scan. A ded- 
icated ring finding algorithm to select these weak rings 
was developed. The vertex and lepton track are very well 
known at this stage of the reconstruction: as single ring 
events, the refined vertex fitter (see section III B I was 
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SK-II Data Total MC Signal MC 

FC, FV, two-ring, 0.1 < Evis < 4 GeV 1023 (100%) 910.8(100%) 26.2 (100%) 
l.KP < 7/6L(m) + 1.1 GeV/c 471 (46.04%) 417.1(45.79%) 18.3(68.8%) 
Angle of 2"^^ cone < 34° 74 (7.23%) 72.73(7.99%) 13.7(52.3%) 
Pattern ID cut 51 (4.99%) 49.64 (5.5%) 12.84 (49.01%) 
-0.75 <V^ < 1.5 GoV^ 12 (1.17%) 18.99 (2.09%) 11.22 (42.81%) 



TABLE VI: Summary of SK-II CCQE selection for two-ring events. In this table the Monte-Carlo has been reweighted according 
to live time, solar activity for SK-II, as well as neutrino oscillations assuming Am^ = 2.5 x 10"'^ eV^ and sin^ 29 = 1. 




FIG. 14: Distribution of analysis variables with cuts for 2-ring CCQE event selection in SK-I. The top two plots are two- 
dimensional L, P distributions for data and MC, shown with the charged pion removal cut; the arrows show the selected 
domain. In the bottom three plots, the full and dashed line show Monte-Carlo expectation from NEUT and NUANCE, with 
oscillation reweighting, but no correction for the absolute flux normalization. The hatched regions show the contribution from 
CCQE events. 



applied, and the e-like or /i-like nature of the track is 
known. The observed charge pattern is projected into a 
50 X 50 grid in (cos 9, 0) spherical coordinates around the 
lepton track. The expected light from the lepton is sub- 
tracted, and the maximum amount of charge in a 25 x 10 
sliding window is found, avoiding directions too close to 
the lepton track. The corresponding average direction is 
used as a starting position in a precise grid fit, where a 
ring-shaped pattern is fitted to the observed light pat- 
tern, varying the direction and opening angle until the 



best match is found. This direction fit yields a new ring 
candidate, which could be a proton or another particle, or 
just a region with higher charge in a bona-fide single ring 
event. 

Therefore, we also calculate the distribution of charges 
Q{9) as a function of opening angle from the candi- 
date direction, subtracting the expected light contribu- 
tion from the lepton. We then calculate the ratio of the 
integrated charge around the maximum of this distribu- 
tion to an average of the charge at low and high values 
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FIG. 15: Distribution of analysis variables with cuts for 2-ring CCQE event selection in SK-II. The top two plots are two- 
dimensional L, P distributions for data and MC, shown with the charged pion removal cut; the arrows show the selected domain. 
In the bottom three plots, the full and dashed line show Monte-Carlo expectation from NEUT and NUANCE, with oscillation 
reweighting, but no correction for the absolute flux normalization. The hatched region shows the contribution from CCQE 
events. 



of 6. This estimator Qdcns, should be large for a true 
Cherenkov ring, since a real ring should have a marked 
peak in Q{9) near the Cherenkov edge. A similar tech- 
nique is also used in standard ring counting to reject fake 
seeds found by the Hough transform. 

Using the lepton's fitted information and the new fit- 
ted direction, we add a lepton and a proton's expected 
charge patterns. The likelihood of this combined pat- 
tern is maximized as already outlined previously, yield- 
ing log(Liopton-i-proton)- Another pattern likelihood cor- 
responding to the lepton alone log(Licpton aiono) is also 
calculated, and the pattern ID estimator 

AL = log [log(Liopton+proton) " log(il cpton alone jj 

is constructed. A high value of this estimator corresponds 
to a proton- like track candidate. 

Our analysis of single-ring events closely follows the 
steps outlined for the two-ring search. 

1. We select single-ring fully-contained events within 
the fiducial volume and 0.1 < Evis < 4 GeV. 



2. Using the output of the proton pattern fitter on 
the new candidate ring, we apply the cuts A < 
P < B X L{m) + C, with A = C =1.1 GeV/c and 
B = 7/6 GeV/c/m. We also require that L < 2 
m. These cuts aim at rejecting rings coming from 
pions as in section |V A[ as well as fake rings which 
are usually fitted with a very low P and L. 

3. We then require logQdeiis > Oi to remove bad ring 
candidates. 

4. In order to improve proton selection we then ap- 
ply the cut AL > 3, where AL is the pattern ID 
estimator defined above. 

5. Finally to reject non-CCQE background as well 
as other non proton rings we ask that > 
—0.75 GeV/c^, using the same definition of V as 
in the previous section. 

The final signal efllciency is about 21% (18% in SK-II), 
for the cases where the proton ring is potentially visi- 
ble. This is lower than in the two-ring case because for 
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single-ring events, most of those protons are just above 
Cherenkov threshold and therefore very hard to find. 
Signal-to-noise after all cuts is about 1.1 to 1. We expect 
to recover 20 (9) CCQE events in SK-I (SK-II), almost 
doubling our statistics. The efficiencies are summarized 
in tables |VII| and |VIII| and the distributions of the vari- 
ables are shown in Figs. [16] and [17] 

In the final sample in SK-I (SK-II), 31 (13) events out 
of the final 39 (18) do have a proton track missed by 
the ring finder that was recovered by this new fitter; 11 
(4) of those are non-CCQE interactions with a visible 
protons (mostly CC pion production) and are an irre- 
ducible background to the CCQE sample. After all cuts 
the single-ring sample contains 90% neutrinos and 10% 
anti-neutrinos. 



C. Kinematic reconstruction of atmospheric 
neutrinos 

Once we have a CCQE enriched sample, kinematic re- 
construction of the incoming neutrino is straightforward, 
assuming that the target neutron is immobile. Writing 
the total momentum as Ptot and the direction as d the 
following equations hold: 



and 



-(Pp+Pi), 



where Pp and Pi are the proton and lepton 3-momenta. 
Since we assume that all selected events are CCQE we use 
E = Ptot for the neutrino energy. Using simulated data 
we have checked the angular resolution on neutrino track 
direction, as well as the energy resolution. The energy 
reconstruction process has a tail toward negative values 
for non-CCQE events reconstructed with this technique 
and is non-Gaussian. We have computed a 68.3% confi- 
dence interval centered on the peak of the distribution. 
The angular and energy resolutions arc summarized in 
table EXl 



1. Zenith angle distributions 

To examine the evidence of zenith angle distortion in 
the enhanced CCQE sample, we first separate it into elec- 
tron and muon neutrino interactions. Each event is cat- 
egorized according to the lepton's e-like or mu-like iden- 
tification. After this characterization, we kinematically 
reconstruct each event. The reconstructed neutrino en- 
ergy spectra are shown in Fig. [18] 

Figure [19] shows the reconstructed neutrino zenith an- 
gle distribution, along with the MC expectation with and 
without oscillations. Clearly, the distributions are consis- 
tent with our the measurement of oscillation parameters 
presented in [5]. 



For the events selected in this analysis, which are 
mostly in the sub-GeV region, the lepton angular cor- 
relation to the incoming neutrino is weak 80° resolu- 
tion). However, with neutrino kinematic reconstruction, 
the angular resolution on the neutrino track reaches 12° 
for and 16° for v^, making the expected zenith an- 
gle distortion in the neutrino direction of mu-like sample 
clearly visible. The observed up-down asymmetry using 
neutrino track reconstruction is — 0.52 ± 0.17(stat) for /x- 
like events and —0.06 ±0.24 for e-like events, compatible 
with Vt neutrino oscillations. For comparison pur- 

poses the asymmetries obtained using the lepton zenith 
angle are -0.059 ± 0.24 for /i-Ukc and ± 0.22 for c-likc 
events. 

The distributions in Fig. [19] can be contrasted with 
those in Fig. 
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In the second set of figures, the lep- 
ton direction rather than the kinematically reconstructed 
neutrino direction is shown for the CCQE enhanced sam- 
ple. As can be seen, there is little sign of the oscillation 
distortion in the lepton zenith angle distributions. The 
correlation between the incoming neutrino and outgoing 
lepton direction has been lost. Most of the events in the 
CCQE enriched sample have a large angle between the 
neutrino and lepton by virtue of the fact that a large 
momentum was given to the proton in the reaction, re- 
sulting in the proton being visible and the lepton having 
a large transverse momentum to balance it. 

Approximately one half of the sample is made up of 
two-ring events where one of the rings has been identi- 
fied as a proton. In order to obtain a high charged current 
purity, in |2] only multi-ring events with a leading muon 
(not an electron) were included. For the same reason, 
a minimum requirement of 600 MeV/c was imposed on 
the leading muon momentum for the sub-GeV multi-ring 
muons. Both of those restrictions have been lifted here. 
In the analysis of |5], even applying the extra require- 
ments used in the mu-like multi-ring sample resulted in 
a CC purity of 54% in the e-like multi-ring sample. The 
use of the proton tag in the CCQE enhanced sample re- 
sults in charged current purity of 88% for e-like events 
and 95% for mu-like events. 

In the combined single ring and two-ring enriched 
CCQE sample the zenith distortion is clearly seen. Ad- 
ditionally, although the statistics is too low to be usable 
(there are only about 20 events in total), according to 
Monte Carlo, the data is expected to show a zenith an- 
gle suppression even for events with the momentum of 
the lepton less than 400 MeV. As shown in 2 , in previ- 
ous analyses, using the lepton direction for these events 
resulted in the complete loss of the initial neutrino direc- 
tion. 

The use of the CCQE sample restores zenith angle 
pointing to a low energy portion of the atmospheric sam- 
ple which previously had no observed zenith angle distor- 
tion, and adds events to the sample with good pointing 
which were not previously considered. This is an impor- 
tant further cross-check of the both the analysis tech- 
nique and the oscillation hypothesis. 
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SK-I Data Total MC Signal MC 

FC, FV, single-ring, 0.1 < Evis < 4 GeV 7224 (100%) 6523 (100%) 94.4 (100%) 

l.K P < 7/6L(m) + 1.1 and L < 2 m 1456 (20.2%) 1488 (22.8%) 44.8 (47.5%) 

Q density cut 235 (3.25%) 213.1 (3.3%) 22.95 (24.32%) 

Pattern ID cut 76 (1.05%) 66.01 (1.01%) 21.04 (22.29%) 

> -0.75 GeV^ 47 (0.65%) 38.69 (0.59%) 20.33 (21.53%) 



TABLE VII: Summary of SK-I CCQE selection for single ring events. In this table the Monte-Carlo has been reweighted 
according to live time, solar activity for SK-I, as well as neutrino oscillations assuming Am^ = 2.5 X 10"^ eV^ and sin 20 = 1. 



SK-II Data Total MC Signal MC 

FC, FV, single-ring, 0.1 < Evis < 4 GeV 3874 (100%) 3463 (100%) 50.71 (100%) 

l.K P < 7/6L(m) + 1.1 and L < 2 m 702 (18.12%) 625.2 (18.1%) 20.45 (40.3%) 

Q density cut 193 (4.98%) 175.2 (5.1%) 12.75 (25.14%) 

Pattern ID cut 29 (0.75%) 27.53(0.80%) 9.58 (18.9%) 

> -0.75 GeV^ 22 (0.57%) 17.6 (0.51%) 9.2 (18.1%) 



TABLE VIII: Summary of SK-II CCQE selection for single ring events. In this table the Monte-Carlo has been reweighted 
according to live time, solar activity for SK-II, as well as neutrino oscillations assuming Am^ = 2.5 x 10~^ eV^ and sin^ 26 — 1. 




FIG. 16: Distribution of analysis variables with cuts for single-ring CCQE event selection in SK-I. The top two plots are 
two-dimensional L,P distributions for data and MC, shown with the charged pion removal cut; the arrow shows the selected 
domain. In the bottom three plots, the full and dashed line show Monte-Carlo expectation from NEUT and NUANCE, with 
oscillation reweighting, but no correction for the absolute flux normalization. The hatched region shows the contribution from 
CCQE events. 



2. L/E distributions uncertainty: 

The neutrino path length L is a simple function of 
the zenith angle however, there are two main sources of 
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FIG. 17: Distribution of analysis variables with cuts for single-ring CCQE events selection in SK-II. The top two plots are 
two-dimensional L,P distributions for data and MC, shown with the charged pion removal cut; the arrow shows the selected 
domain. In the bottom three plots, the full and dashed line show Monte-Carlo expectation from NEUT and NUANCE, with 
oscillation reweighting, but no correction for the absolute flux normalization. The hatched region shows the contribution from 
CCQE events. 



two-ring events single-ring events 
CCQE TOTAL CCQE TOTAL 



SK-I 


Angular resolution (°) 
Energy resolution (%) 


8.0° 
+ 10.7 

- 4.4 


12.3° 

+ 10.6 
- 17.2 


8.1° 
+ 9.9 

- 6.2 


14.9° 

+ 13.6 
- 16.7 


SK-II 


Angular resolution (°) 
Energy resolution (%) 


8.6° 

+ 10.3 
- 5.9 


12.0° 
+ 10.8 
- 17.9 


8.1° 

+ 9.6 

- 5.8 


16.2° 

+ 10.8 
- 22.0 



TABLE IX: Energy and angular resolution of the kinematic reconstruction shown in the previous paragraph. Since the energy 
resolution distributions are asymmetric we report 68.3% centered confidence intervals. 



• Near the horizon, a small error on the zenith angle 
leads to a large error in L since dL / d cos 6 zenith is 
large. 

• There is an additional uncertainty in the produc- 
tion height of atmospheric neutrinos. While it is 
negligible compared to the overall travel length for 
neutrino produced on the far side of the Earth, it 
is comparable to the total L for neutrinos produced 
directly above Super-Kamiokande. The model from 



Honda et a?. [9] is used to estimate the production 
height. For a given zenith angle, L is the average 
path-length over 20 samplings from the production 
height distribution. 

In this analysis, we use the reconstructed neutrino 
zenith angle to estimate L; i? is the kinematically re- 
constructed V energy defined in section |V C| This has 
not been attempted so far in Super-Kamiokande; the 
method reported in |18j used an energy estimator based 
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reconstructed v energy (MeV) reconstructed v energy (MeV) 

FIG. 18: Reconstructed energy spectra of the CCQE enriched sample, e-like events are on the left and mu-like events on 
the right. The hatched region show the CCQE fraction from NEUT. The solid histograms correspond to our non-oscillated 
Monte-Carlo simulation, while the dashed histograms show the prediction with oscillations assuming Am^ — 2.5 x 10~^ eV'^ 
and sin^ 26 ^ I. 
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FIG. 19; Reconstructed neutrino zenith angle of the CCQE enriched sample, e-like events are on the left and mu-like events 
on the right. The box histograms correspond to our non-oscillated Monte-Carlo simulation, while the dashed histograms show 
the prediction with oscillations assuming Am^ = 2.5 x 10~^ eV^ and sin^ 29 = 1. 



30r 



Tagged CCQE e-lil<e 
lepton direction 



g20h 



10: 



30 



g20 
0) 



<D 



'10 



Tagged CCQE mu-lil<e 
lepton direction 



-1 -0.5 0.5 1 -1 -0.5 0.5 1 

lepton cose^^„„, lepton cose^^^„, 

FIG. 20: The lepton zenith angle of the CCQE enriched sample, e-like events are on the left and mu-like events on the right. 
The box histograms correspond to our non-oscillated Monte-Carlo simulation, while the dashed histograms show the prediction 
with oscillations assuming Am^ — 2.5 x 10~^ eV^ and sin^ 26 — 1. 
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on the total energy of all charged particles in the ID, and 
the neutrino zenith angle was obtained using the lepton 
zenith angle. Figure 21 is the distribution of L/E for 
events with a primary /i-likc ring in the CCQE sample. 
SK-I and SK-II and single- and two-ring samples were 
combined together to increase statistics. The agreement 



with the oscillated distribution is very good. 

These distributions can be used in — s- Vt oscilla- 
tion searches. Binning the e-like and ^-like data and MC 
log]^g(i/i?) distributions (5 equally spaced bins from 
to 5 for both samples), we use the following Poisson like- 
lihood ratio [T9l: 



= ^ ( {Muc - Ndata) + Ndata log ^^^'^^ ' 



i=2 



MC 



r 



J=2 



aj is the estimated uncertainty on each Sj parameter; 
these variables describe the sources of systematic errors 
and are adjusted during the fit. Nmc is the weighted sum 
of Monte-Carlo events, taking into account livetime and 



solar activity. Oscillation probabilities and systematic 
errors are taken into account using a technique similar to 
that of 20J. The weight of each Monte-Carlo event is : 



w = F{l + ei){Ej2Ge\Y^P{sin^ 29, Am^, (1 + e3)L/E^){l ± e4/2)(l ± S5/2), 

I 



where is a numerical factor accounting for livetime 
normalization and solar activity, E^, is the true neutrino 
energy, L the neutrino path length, and the function P 
is the oscillation probability. For each event, 20 values 
of L are obtained by sampling the production height dis- 
tribution; P is the arithmetic average of the oscillation 
probabilities obtained for each of these lengths. The ej 
parameters have the following meanings : 

• El represents the uncertainty on the absolute nor- 
malization of the sample. This includes flux un- 
certainties (8% in this energy range [9]), and also 
the systematic uncertainty on our CCQE event se- 
lection, which we estimate to be 14.4% based on 
Monte-Carlo studies. Since the combined error is 
large, it is a free parameter in the flt. 

• £2 is the uncertainty on the neutrino flux spectral 
index, with <J2 = 0.05 as in PU] . 

• £3 is a systematic error on the Monte-Carlo L/E 
ratio, coming from uncertainties in the production 
height, and absolute neutrino energy scale. Follow- 
ing [5] its width is set to 10%. 

• £4 describes the systematic error on e/ fi separation 
and uncertainties on the flavor content of the at- 
mospheric flux: in the definition of w the sign is 
negative for /i-like and positive for e-like, leading 
to complete anti-correlation between the and i>e 
samples. Since SK-I and SK-II are combined, we 
have used a conservative estimate of 15% for this 



• £5 accounts for the uncertainty in the non-CCQE 
background selection efficiency, as well as for the 
CCQE/non-CCQE cross-section ratio; the sign pre- 
ceding £5 in the definition of w is positive for CCQE 
and negative for non-CCQE events. We have as- 
sumed this ratio to be known to 10% in our simu- 
lation. 

• Finally one last parameter is used to shift the 
reconstructed value of logiQ{L/E) to account for 
biases in the kinematic reconstruction technique. 
Combining energy and direction reconstruction er- 
rors, we conservatively estimate the width of this 
parameter to be 10%. 

We have used a 201 x 201 grid in log;^o Am^ and sin^ 29; 
at each point on the grid the ej were fitted using MI- 
NUIT. The best fit point is found in the unphysical 
region at (Am^ = 1.3 lO"'^ eV^ sin^ 26* = 1.3), with 
Xmin — 6.70/8 dof. The best fit in the physical re- 
gion occurs at (Arn^ = 3.410"^ eV^, sin^ 26* = 1.0) at 

shows the 68.3%, 90% 
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xLn^phys = 7-98/8 dof. Figure 
and 99% allowed regions ; these contours are located at 
xliin.phys + 2-86, 5.32, and 10.02 respectively, using the 
same method as in [2 . With the low statistics of this 
sample, the allowed region is large. They are consistent 
with our previous, more precise results: the best fit point 
reported in our previous publications is accepted at bet- 
ter than 1(7. The no-oscillation point is rejected at 3(7 
{xioosc-xlitn,phys = 12.95). The systematic parameters 
at the best fit are summarized in table [X] The kinemati- 
cally reconstructed CCQE sample used here has roughly 
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FIG. 21: Distribution of L/E for the combined CCQE enriched /x-like (left) and e-like (right) samples. The continuous lines 
show the expectation from NEUT without v oscillations, while the dashed lines shows the expectation with Vr oscillations 

with Am^ = 2.5 x 10'^ eV^ and sin^ 261 = 1 



half the statistics of the Soudan-2 tracking calorimeter 
experiment for which a similar analysis was performed 

m- 

It must be noted that the CCQE selection method se- 
lects a quasi-pure sample of neutrinos (~ 90%). This 
could potentially have sensitivity to the mass hierarchy 
and other CP-odd effects. 



VI. CONCLUSION 

We have developed a technique for identifying pro- 
ton tracks in the Super-Kamiokande detector, relying 
on their distinctive Cherenkov ring pattern. The pro- 
ton momentum threshold for Cherenkov radiation is 1.07 
GeV/c, and above about 2 GeV/c they produce visible 
secondary particles through hadronic interactions in the 
water, which makes identification outside of this proton 
momentum range very difficult. Because of this we find 
that protons can be identified efficiently in the momen- 
tum range between 1.25 GeV/c and 1.7 GeV/c. The 
majority of atmospheric neutrino interactions produces 
protons with momenta below the Cherenkov threshold, 
which results in low statistics. Using the full SK-I and 
SK-II data sets (2285.1 days), we have applied this tool 
to identify neutral current elastic events. The com- 
bined data is compatible with the observation of protons 
(X^ = 9.3 for 6 bins, P = 15.7%). 

In a separate analysis, we have selected a sample of 
CCQE events, reconstructing both the lepton and the 
proton, and thereby the incoming neutrino kinematic pa- 
rameters. This is the first time that kinematic recon- 
struction of the neutrino track is reported in a water 
Cherenkov detector detector using atmospheric neutri- 
nos. This sample has 91.7 ± 3% neutrino (as opposed to 
anti-neutrino) content, potentially allowing studies where 
separating neutrinos and anti-neutrinos is crucial, e.g. 
resolving the neutrino mass hierarchy. The technique 



reported in this paper has an angular resolution of 12° 
for and 16° for on the neutrino track, allowing 
the observation of a clear zenith-angle dependant asym- 
metry in the neutrinos themselves. The observed up- 
down asymmetry using neutrino track reconstruction is 
-0.52 ±0.17 (stat) for ^-like events and -0.06 ±0.24 for 
e-like events, compatible with our previous results with 
other higher energy samples and — > neutrino oscil- 
lations. 

Although the CCQE selection yields low statistics in 
Super-Kamiokande and is not competitive with previ- 
ously published measurements of the oscillation param- 
eters, we have performed a fit of the L/E distribution 
including relevant systematic effects. With this sam- 
ple alone, the no-oscillation hypothesis is excluded at 
3 standard deviations. Table IXII summarizes the new 
data samples selected in this paper. We expect that 
this new technique will be of particular interest to fu- 
ture projects involving water Cherenkov detectors even 
larger than Super-Kamiokande, where the statistics will 
be much higher. 
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FIG. 22: Allowed regions at 68.3% (dotted), 90% (dashed) and 99% (continuous) CL for — > v-r oscillations from the CCQE 
kinematically reconstructed samples. The star shows the positions of the best fit from [2], while the triangles show the position 
of the best fit for this sample in the physical region. These allowed regions are consistent with previously published results and 
exclude the no-oscillation hypothesis at greater than 3 sigma. 



parameter meaning uncertainty value at best fit 



El absolute normalization (free) 6.5% 

£2 spectral index 0.05 -0.0006 

£3 Error on true L/E 10% -1.9% 

£4 e/^ ratio 15% -2.1% 

£5 CCQE/non-CCQE ratio 10% 0.2% 

£6 shift in reconstructed L/E 10% -5.3% 



TABLE X: Systematic parameters and their best fit values. 



Event clciss SK-I data SK-II data 

NC elastic (expected NC elastic fraction) 27 (64.7%) 11 (55.6%) 

CCQE e-hke (expected CCQE fraction) 31 (53.0%) 16 (51.4%) 

CCQE /x-like (expected CCQE fraction) 60 (62.4%) 18 (61.3%) 



TABLE XL Summary of the observed data samples : single proton, CCQE e-like and CCQE /i-like, for SK-I, SK-II. 
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